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ABSTRACT 



We have mapped the central lO'xlO' of the dense core TMC-IC at 450, 850 
and 1200 /xm using SCUBA on the James Clerk Maxwell Telescope and MAMBO 
on the IRAM 30m telescope. We show that although one can, in principle, use 
images at these wavelengths to map the emissivity spectral index, temperature 
and column density independently, noise and calibration errors would have to 
be less than ~2% to accurately derive these three quantities from a set of three 
emission maps. Because our data are not this free of errors, we use our emission 
maps to fit the dust temperature and column density assuming a constant value 
of the emissivity spectral index and explore the effects of noise on the derived 
physical parameters. We find that the derived extinction values for TMC-IC are 
large for a starless core (~80 mag Ay), and the derived temperatures are low 
(~6 K) in the densest regions of the core, using our derived value oi (3 = 1.8. 

Subject headings: stars: formation — dust, extinction — submillimeter 



1. Introduction 



Efforts to determine the mass and temperature of starless cores from sub-millimeter 
and millimeter observations are hampered by uncertainties in the emission properties of 
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the dust grains, such as the emissivity spectral index. Although in principle it should be 
possible to calculate the column density of dust, the emissivity spectral index of the dust, 
and the dust temperature from observations at three or more wavelengths, in practice this 
has never b een done for a starless cor e . A similar analysis has been done for circumste llar 
disks, e.g. (IBeckwith fc Sargentlll99ll : [Mannings fc Emerson! 1 1994 : iMathieu et al.lll995l ) in 
which temperature gradients, disk masses and spectral indices are calculated. In this paper 
we explore the levels of uncertainty in the derived dust temperature (T^), emissivity spectral 
index (/?) and column density (Nh) resulting from datasets of either three or four noisy 
emission maps at different wavelengths. We then apply this analysis to the starless core 
TMC-IC. 



TMC- 



C is a starless core in the Taurus molecular cloud, at an approximate distance 



of 140 pc (IKenyon et al.l 119941 ). It was shown that TMC-IC is a coherent core, meaning 



that its velocity di spersion is roughly constant, at slightly more than the sound speed, over 
a radius of 0.1 pc (IBarranco fc Goodmanlll998l : iGoodman et al.lll998l). The velocity field o f 



TMC-IC shows evidence of solid body rotation, at 0.3 km s ^ pc ^ ( Goodman et al. Il993 ). 



and the N2H"'"(l-0) spectrum reveals the signature of sub-sonic infall (ISchnee &: Goodman 
20051 ). The mass derived from 450 and 850 /xm maps alone (13 Mq) is se veral times the virial 
mass , and the density profile is similar to that of a Bonnor-Ebert sphere (ISchnee &: Goodman 
2005h . 



Here we use data taken with SCUBA (at 450 and 850 /xm) and MAMBO (at 1200 /jm) 
to make maps of the dust column density and temperature, and to estimate a constant value 
for the emissivity spectral index of TMC-IC. Although the high signal to noise at 850 and 
1200 /im make this set of maps one of the best yet available for a starless core, we show 
that the noise is still too high to reliably map variations in the emissivity spectral index of 
TMC-IC. 



2. Observations 



2.1. SCUBA 



We observed a lO'xlO' region around TMC-IC using SCUBA flHoUand et al.lll999[ ) on 
the JCMT. Our maps, especially at 450 fim, benefitted from exceptionally stable grade 1 
weather. We us ed the standard scan-rnapping mode, record ing the 850 and 450 /xm data 
simultaneously (jPierce-Price et al.ll2000l : iBianchi et al.l 120001 ). Three chop throw lengths of 
30", 44", and 68" were used in both the right ascension and declination directions. The 
JCMT has FWHM beams of 7.5" at 450 /xm and 14" at 850 /im, which subtend diameters of 
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0.005 and 0.01 pc, respectively, at the distance of Taurus. Pointing during the observations 
was typically good t o 3^^ or better. The data reduction for the SCUBA data is described by 
Schnee &: Goodmaru (120051 ). The absolute flux calibration is uncertain at levels of ~4% at 
850 iJ,m and ~12% at 450 fim. The rms noise in the 850 /im map is 8.6 mJy/beam, and 13 
mJy/beam in the 450 /im map, measured in regions with no significant emission. 

The Emerson2 technique used to reconstruct SCUBA sca n maps from its comp onent 
chop throws introduces false structure that can be removed (jjohnstone et al.ll2000l ). To 
remove this structure, we convolved the SCUBA images (after masking out pixels with 
\S\ > 5a) with a Gaussian of FWHM twice t he size of the largest ch op throw, and subtracted 
this from the original image, as explained in lReid fc Wilson! (120051 ). The resulting image has 
fluxes nearly identical to the original in regions of high signal to noise, but has fewer bowls 
of negative emission and other artifacts introduced by chopping and image reconstruction. 
We only use data from the high signal to noise central region of the maps to estimate the 
dust properties. We account for the SCUBA error beams by convolving the 850 and 1200 
/im maps with the 450 fim PSF, and convolving the 450 and 1200 /im ma ps with the 850 
um P SF before regridding to a common resolution, as explained in detail in iReid &: Wilson 



The original bowls of negative flux in our 450 fim map on either side of region 'a' in 
Figure[7]have peak values around —150 mJy, and average values around —75 mJy. We believe 
that these structures are artifacts of image reconstruction. It is possible that structures of 
similar magnitude (either positive or negative) might be affecting data elsewhere in the map. 



2.2. MAMBO 

We observed the 1.2 mm continuum emission in November 8, 2002, October 23 , 2003, 



and November 2, 2003 using the 117-channel MAMBO-2 array (iKreysa et al.lll999l ) at the 
IRAM 30-meter telescope on Pico Veleta (Spain). The FWHM beam size on the sky was 
10'.'7. The source was mapped on-the-fiy, with the telescope subrefiector chopping in azimuth 
by 60" to 70" at a rate of 2 Hz; the total on-target observing time was about 6 hours. The 
line-of-sight optical depth varied between 0.1 and 0.5. The data were reconstructed using the 
EKH algorithm in an iterative way that properly reproduces large-scale emission (Kauffmann 
et al., in prep.). The rms noise in the 1200 /im map is 3 mJy/beam, measured in regions 
with no significant emission. The flux calibration uncertainty is approximately 10%, which 
is derived from the rms of calibrator observations across pool observing sessions and the 
uncertainty in the intrinsic calibrator fluxes. 
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3. Solving for physical parameters 

In principal, one can use three measured quantities (e.g. the three flux density maps 
at 450, 850 and 1200 /xm) to solve for three unknowns (e.g. maps of the dust temperature, 
emissivity spectral index, and column density). Our maps originally are in units of flux 
density per beam, with a different beam size for each map. To derive meaningful physical 
quantities we smooth and rebin our maps to a common resolution of the largest beam, which 
in our case is 14". As a result, the quanity that we work with is flux density per 14" pixel, 
which is what we present in our maps of TMC-IC. 

The flux density per beam in each map is given by: 

= nB^{Td)K^lj,mHNH2, (1) 

where 

r 

R (TA = - 

exp(/iz//A;Tc 



. , 2hu^ 1 
BAT,) = . (2) 



and 

\ /3 



= "^^4230GIfe j ■ 

In Equation [H is the flux density per 14" pixel; Q is the solid angle of the beam; ByiTd) 
is the blackbody emission from the dust at temperature T^; K^.^n = 0.0 09 cm^ g~^ is the 



emissivity of the dust grains at 230 GHz (jOssenkopf fc Henningl Il994l ): mij is the mass 



of the hydrogen atom; ^ = 2.8 is the mean molecular weight of interstellar material in a 
molecular cloud per hydrogen molecule; is the column density of hydrogen molecules and 
a gas-to-dust ratio of 100 is assumed. It should be noted that ^230 is uncertain by a factor 
of ~2, and that we assumed that K230 = 0.005 and = 2.33 (which is the mean molecular 
w eight per free particle for an abundance ratio of N{H)/N{He) = 10 and negligible metals) 
in ISchnee fc GoodmanI (120051 ). 

The ratio of two fluxes, because of the common beam size, can be simply expressed as: 

Su2 V^2/ \exp[hiyi/kTd] - 1 J 

The dust temperature can be found independently of the dust emissivity spectral index 
by taking the difference between the ratio of fluxes, if we assume that each line of sight 
through the core can be characterized by a single temperature and emissivity spectral index: 



, ,^450 A, /1200/im\ 
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1 / 'S'sso \ , f 850/im 
log 7; log 



5'i2oo / v450/xm 
, /exp [AT/850uml - 1\ , /1200um\ 

log rx lArn -^ 7 log 



exp [Ar/450/im] — 1/ \ 850;um y 
/exp [AT/1200/im] - 1\ /850/im\ 
V exp [AT/850;um] - 1 J ) ' 

where At = hc/kTa. 

Once the dust temperature is determined, the emissivity spectral index can be calculated 



by: 



exp [AT/850/.m] - 1 \ / / 1200^\ _ ^ 
\S12m exp [AT/1200/im] - 1/ / V 850/xm / ^ ^ 



The column density of dust can be derived from the flux at a single wavelength (e.g. at 
1200 yLim), the temperature of the dust and the emissivity spectral index of the dust using 
Equation [H The equivalent visual extinction can be calculated from the column density 
using: 

Ay = NuRv^^^-^ (7) 

where Nh = 2 x Nh2i Nh/E{B — V) = 5.8 x lO^""^ cm~^ mag""*^ is the conversion between 
column density of hydrogen nuclei (for our assumed gas to dust ratio) and the selective 
absorption, and Ry = Ay/E{B — = 3.1 is the total to selective extinction ratio for the 



low-density lines of sight similar to those for which Nh/E{B~V) has been measured (IMathis 



I990I : iBohlin et al.lll978l ). Although we use a constant value of Ry, we recognize that this 



value is uncer tain to within a factor of 2 between regions of high and low column density 



( lMathislll990l ). and the relation between extinction and column density may be different for 



dense cores like TMC-IC. 



4. Error Analysis 

In order to understand how the noise and reconstruction artifacts in our emission maps 
will affect the accuracy of our derived dust temperature, emissivity spectral index and column 
density we have run a variety of Monte Carlo simulations. We compare the effects of using 
maps at three particular wavelengths to solve for all three parameters to using the three 
measurements to fit for two parameters while assuming a fixed value for the third. We also 
show the improvements brought about by using a fourth wavelength and fitting for all three 
physical parameters. 
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4.1. Illustrative Examples 

Here we discuss the uncertainties to be expected from solving for T^, (3 and Nh from 
three emission maps that have noise or whose calibration is uncertain. To illustrate the 
method, Figure [1] shows the modified blackbody spectrum from dust with Ay = 50, = 10 
K and (3 = 1.5, observed with a 14" beam. The black dotted line shows the true emission 
spectrum, with crosses at 450, 850 or 1200 fim, and the blue/red crosses show the flux 
overestimated/underestimated at one wavelength by 20%. The blue and red curves are the 
fitted spectra that pass through the one blue/red cross and the two black crosses. It is clear 
that a 20% error in one measurement creates errors in all three derived parameters that are 
much larger than 20%. The derived values of T^, (3 and are labeled in Figure [TJ For 
convenience, we show the derived Nh in units of Ay- 

As can be seen from Equation [5], noise in any of the three flux maps first translates into 
uncertainty in the derived temperature. This incorrect value of T^, along with any noise 
in the 850 and 1200 fim maps then results in uncertainties in the derived values of (3 and 
Nh- From Equation [6l one can show that overestimating the temperature will result in an 
underestimate of /3, and underestimating the temperature will result in an overestimate of f3. 
The errors in the derived parameters from incorrectly measuring the flux at one wavelength, 
while correctly measuring the flux at the other two wavelengths is shown in Figure [H The 
anti-correlation between and (3 is clearly seen. Also apparent is that for dust with "core- 
like" values of (3, Td and Nh-, errors in 450 /im flux result in smaller errors in the derived 
physical parameters than errors at 850 and 1200 /im, which is convenient because the 450 
/im maps often suffer from higher levels of noise than maps at 850 and 1200 /xm. 

4.2. Deriving the Physical Parameters 

In order to determine the effect of similar levels of noise in all three wavelengths on 
the three derived physical parameters, we determine the flux at 450, 850 and 1200 /im from 
dust at Td = 10, 15 and 20 K and (3 = 1.0, 1.5 and 2.0. The fluxes are then modified by 
a multiplicative factor / = 1.0 + 5 where 6 is randomly chosen from a normal distribution 
of mean zero and standard deviation a, and each flux is modified by a different /. This is 
repeated 10,000 times for each value of a between and 0.2. We show the effect of noise 
in all three wavelengths on the derived column density, temperature and emissivity spectral 
index in Figure [3l At a signal to noise of 20 (5% error), the expected uncertainties in the 
derived Nh, Td and P are approximately 50%, 80% and 40%, respectively, for 15 K dust 
with f3 = 1.5. The median values of the derived parameters stay close to the input values at 
every value of Td and /3 tested. 
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The analysis presented here deals with the special case that the signal to noise is wave- 
length independent. In general, this will not be the case, and one might expect that for a 
given amount of observing time the signal to noise will be worse at 450 fim than at 850 or 
1200 fim due to atmospheric effects. In addition, the relative signal to noise between maps 
will change from position to position, due to the gradients in the dust temperature, column 
density and emissivity spectral index. For instance, at the position of the column density 
peak in TMC-IC, the signal to noise ratios at 450, 850 and 1200 /im are 5.3, 27 and 33, 
respectively, when including both the random noise and the artifacts in the 450 fim image. 

Although even just 5% errors in the measured fluxes at three wavelengths make accurate 
determinations of three physical parameters impossible, adding a fourth wavelength (for 
instance, at 350 /xm or 2.7 mm) drastically reduces the effects of noise, as can be seen by 
comparing Figure [3] with Figure |H To properly constrain the dust parameters, at least one 
of the observations should not be on the Rayleigh- Jeans portion of the emission spectrum, 
or else and Nh will be degenerate in the fit. For dust at Td = 10, 15 and 20 K and 
f3 = 1.0, 1.5 and 2.0 we calculate the flux at 350, 450, 850 and 1200 /im. As before, the 
fluxes are then modified by a multiplicative factor / = 1.0 + 6, where 6 is again a variable 
randomly chosen from a normal distribution with mean zero and standard deviation a. This 
is repeated 10,000 times for each value of a between and 0.5. We allow a to be larger 
than in the previous Monte Carlo simulation because the effects of noise are smaller in this 
case. For each set of four fluxes, the column density, temperature and emissivity spectral 
index are fit and the results are shown in Figure |H At a signal to noise of 20 (5% error), the 
expected errors in the column density, temperature and emissivity spectral index are all on 
the order of ~1%. As the signal to noise gets lower, the median derived temperature and 
column density decrease, and this effect is larger for warmer cores. 



4.3. Fixing One Parameter 

In the case of a starless core observed at three wavelengths, one can hold one parameter 
fixed, such as assuming a constant = 15 K or f3 = 1.5, and use the three flux measure- 
ments to fit the remaining two parameters. Figure [5] shows the result of using noisy flux 
measurements to fit the column density and either the dust temperature or the emissivity 
spectral index of a core with = 15 K and f3 = 1.5. We show the results for correctly 
assuming that /? = 1.5 and erroneously assuming that /3 = 2.0 or 1.0, as well as for correctly 
assuming that = 15 K and erroneously assuming that = 10 or 20 K. 

We see that by choosing the correct value of (3 and fitting the Td and Nh, noise in the 
emission maps on the level of 10% result in ~1% uncertainties in the derived temperature 
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and column density. Choosing a value of /5 = 1.0, when the proper value is /3 = 1.5 results 
in temperatures that are on average 3 K too high, with a spread of ~3 K, and column 
densities that are high by 8%. Overestimating the emissivity spectral index, by assuming 
that 13 = 2.0, results in temperatures that are on average 2 K too low and column densities 
that are 5% low. 

Figure [5] shows that by choosing the and fitting for the (3 and Nh, the true values of 
/3 and column density are recovered, on average, with an absolute spread of 0.20 in f3 and 
a relative spread of 13% in column when the signal to noise ratio is 10. Underestimating 
the temperature by assuming = 10 K when the proper value is = 15 K results in 
an overestimate of the emissivity spectral index by 0.6, on average, and an overestimate of 
the column density by 70%. Overestimating the temperature by assuming that T^; = 20 K 
results in an underestimate of (3 by 0.3, on average, and an underestimate of Nh by a factor 
of 30%. 



4.4. A Model Core 

We apply the error analysis presented in Section H] to our data on the starless core 
TMC-IC in order to derive new sci ence beyond what was possible in our earlier paper on 



this core (ISchnee &: Goodman! l2005l ). The random noise in the emission maps, as measured 
in regions with faint emission, is found to be 13 mJy/beam, 9 mJy/beam and 3 mJy/beam, 
which corresponds to S/N values of 21, 14 and 15 at 450, 850 and 1200 fim at the (0,0) 
position of the emission maps. Calibration uncertainties are ~12%, 4% and 10% at 450, 850 
and 1200 /im, respectively, and the image reconstruction artifacts in the 450 map peak 
at 150 mJy/beam. 

To determine how the random noise and reconstruction artifacts (which are spatially 
correlated and therefore not truly random) in the observed emission maps affect the derived 
parameters, we construct synthetic emission maps of a starless core like TMC-IC. The model 
core has a temperature and column density profile equal to the one derived from the two 
dimensional temperature and column density profile of TMC-IC, derived from fitting the 
450, 850 and 1200 /im data and assuming a constant f3 = 1.5. The model starless core is 
made of cylindrical shells seen face on with a central temperature of ~6 K, rising to ~12 K 
at the edge. The central column density corresponds to ~80 magnitudes of visual extinction, 
falling to an Ay of ~20 at the edge. Using the equations in Section [3l we derive the resultant 
emission maps at 450, 850 and 1200 fim, add in Gaussian noise of the same magnitude as the 
random noise and reconstruction artifacts described above, and from these emission maps 
derive maps of the column density, temperature and emissivity spectral index. The resultant 
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maps are shown in Figure [6l 

The dust emission in our model core (shown in Figure [H]) at 450 and 850 /xm is not a good 
tracer the dust column density, though the 1200 fim emission map does resemble the column 
density distribution. Spatial gradients in the dust properties, such as the temperature and 
emissivity spectral index, need to be taken into account when attempting to find the "peak" 
of the dust distribution, even when using the relatively longer wavelength 1200 fim emission 
map as a proxy for column density. 

When we solve for all three physical parameters, we find that the derived (3 of our 
model cloud has a median value of 1.5 (which is the input (3 everywhere in the model) with a 
standard deviation of 0.7. Given the close correspondence between the input value of (3 and 
the median value derived for it, we see that a constant value for the emissivity spectral index 
can be estimated in this manner. Using this constant value for j3 everywhere, we can then 
use our three fiux maps to fit the temperature and column density. The resultant maps of 
our model constructed in this way are shown in Figure [61 Although using the median value 
in the /3 map to derive a constant value for the emissivity spectral index reduces the impact 
of statistical uncertainty due to random noise, systematic shifts, such as those created from 
calibration uncertainties, are not removed. 



5. Dust Emission in TMC-IC 
5.1. Morphology 

The observed 450 /im emission map of TMC-IC is qualitatively different from the 1200 
/xm map. The 450 fim map shows a condensation at (50, —150) (position "a" in Figure 
[7]) which is much fainter at 850 /xm and nearly absent at 1200 fim. The condensation at 
(—50, 150) (position "c" in Figure [7]) is prominent at all wavelengths, while the column 
density peak at (0,0) (position "b" in Figure [7j) is prominent at 850 and 1200 /xm, but not 
apparent at 450 /im. 

An emission peak in the longer wavelength maps, but not prominent at 450 /im, can be 
explained by a cold temperature, as seen at the (0, 0) position. The 450 /xm emission peak 
not seen at 1200 /im can be explained by the dust in that region having a steep spectral 
index (/? > 2). 
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5.2. Derived Parameters 



Based on our analysis in Section I4.H calibration uncertainties, reconstruction artifacts 
and noise in our TMC-IC emission maps prevent us from making accurate maps of dust 
temperature, emissivity spectral index and column density simultaneously, even though this 
may well be the highest S/N set of such maps of a starless core to date. Figure [7] shows the 
results of an attempt to do so, and as expected the range of temperatures that we see is 
quite broad (T^ < 5 K and Td > 25 K) for a starless core, with a similarly large spread in 
emissivity spectral index (/? < and f3 > 2.5) and column density. Furthermore, the errors 
and uncertainties in our emission maps create a spurious anti-correlation between the dust 
temperature and emissivity spectral index, which is also seen in our attempt to derive T^, 
j3 and Nh from our simple model of a cylindrically symmetric starless core, described in 
Section 14. 4[ Because the noise and errors in our observed emission maps of TMC-IC will 
drive an anti-correlation between the derived dust temperature and emissivity spectral index 
even when such a trend does not exist, our results are consistent with a constant value of 
the emissivity spectral index. However, we cannot rule out a real anti- correlation such as 
that observed at 200, 260, 360 and 5 80 fj,m in the Orion molecular cloud (IDupac et al.ll200ll ) 
and the Ml 7 star-forming complex (IDupac et al.ll2002l ). 



Following the method described in Section we take the median value of the emissivity 
spectral index map (Figure [7]) and derive a value and uncertainty of /5 = 1.8 ± 0.5 for the 
portion of TMC-IC located within the white contour of the 850 fim image in Figure [71 
Unless otherwise state, in our subsequent analysis of TMC-IC we use a constant value of 
/3 = 1.8 everywhere in the core. We derive (3 and the la uncertainty by running 1000 
realizations of a Monte Carlo simulation of the observed flux maps in TMC-IC modified by 
the calibration uncertainties. Using /5 = 1.8, we construct dust temperature and column 
density maps (Figure [8]) from a fit to the 450, 850 and 1200 /im images. The column density 
in Figure [8] peaks around the maximum of the 850 and 1200 //m emission. The implied visual 
extinction is quite high, rising above 80 magnitudes in the densest regions. As expected, the 
regions with the high est column density are also the regions with the lowest dust temperature 
(IZucconi et al.ll200ll ). By using a constant value of the emissivity spectral index {j3 = 1.8), 
the dust temperature that we derive is nowhere significantly higher than 16 K nor lower than 
6 K. 

Our derived value of jS = 1.8 ±0.5 is somewhat h igher than that measured for amorphous 
carbon grains (/? = 1.2) bv iMennella et al.Nl998l) and is within the range (~ 1.2 — 2.5) 
measured for silicate grains by lAgladze et al.l (119961 ) . Dust in int erstellar disks are generally 



observed to have values of 3 < 1, but with con siderable spre a d (IBeckwith fc Sargentlll991 



Mannings fc Emerson 1994). Observations by Stepnik et al. ( 20031 ) have shown that P 
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1.9 ± 0.2 for a dense filament in Taurus, which agrees very well with our estimate of TMC- 
IC. Graphite and silicate dust grains in the ISM are often assumed to have (3 = 2 (e.g. 



Draine fc Led 1 19841 ). 



5.3. Comparison With Previous Results 



Two-dimensional temperature and column density maps of TMC-IC, along with the de- 
pr ojected three-dimens i onal t emperature and density profiles, have previously been reported 



m 



Schnee fc Goodman! (120051 ) using only SCUBA 450 and 850 /im data. With the addition 



of the MAMBO 1200 /xm map, we are better able to constrain the temperature and density 
and estimate the emissivity spectral index. The emiss ivity spectral index that we use here 
(/3 = 1.8) is higher than in ISchnee fc GoodmanI (120051 ) (/? = 1.5), because in this paper we 
are able to derive /3, and in our ealier work we had to assume a valu e. The value that we 
use in this paper for /t o^n is taken from lOssenkopf fc Hennind (119941 ). which is larger than 
the number we used in lSchnee fc GoodmanI (120051 ). As a result using larger (3 and K230, the 
temperatures that we derive are somewhat lower and the densities are also lower. We choose 
here to use the dust opacity appropriate for dust grains with thin ice rnantles , evolved in a 
dense (10^ cm~^) region for 10^ years derived in lOssenkopf &: Hennind (119941 ) because this 
is the consensus value settled u pon by the Spitzer Legacy Project, "From Molecular Cloud 
Cores to Planet Forming Disks" (lEvans et al.ll2003l ). which will make comparisons with other 
c2d cores easier in the future. The mass that we derive for TMC-IC, within a radius of 0.06 
pc fro m the column density peak is 6 Mq, as compared with 13 Mq in 



Schnee fc Goodman 



( I2OO5I ). However, even this new dust-derived mass is higher than the virial mass derived 
from TMC-IC N2H+ observations. 



Following the method used in ISchnee &: GoodmanI (120051 ) . we create deprojected three 
dimensional temperature and density profiles. We assume that the inner 0.07 pc of TMC-IC 
can be approximated as a set of nested spherical shells, each with a constant temperature 
and density. The derived temperature and density profiles are shown in Figure [HI along with 
the profiles calculated usin g only the 450 and 850 /im data and assuming that P = 1.5, as in 
Schnee fc GoodmanI (l2005r). The dust temperatu re profile that we derive is not much changed 
from that derived in ISchnee fc GoodmanI (120051 ) , and is thus still consistent with the dust 
temperature profiles predict ed for extern a lly he at ed starless cores with Bonnor-Ebert density 



distrib utions calculated by lEvans et al.l (1200 ll ): iGoncalves et al.l (120041 ): IStamatellos et al. 
(I2OO4J ). A Bonnor-Eber t profile is a good fit to the TMC-IC density profile at radii greater 



than 0.004 pc, but as in lSchnee Sz GoodmanI (l2005l ). the density we derive for the innermost 
point is significantly higher than predicted by a Bonnor-Ebert model. The density profile. 
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shown in Figure [9|, is consistent with a broken powerlaw, with n{r) oc r~°'^ inside 0.035 
pc and n{r) oc r~^'° outside 0.035 pc. This is considerably flatter than the density profile 
derived using only the 450 and 850 fim data and /5 = 1.5, which have powerlaw exponents 
of —0.8 and —3.1 inside and o utside the break radius, respectively (also shown in Figure [H]). 
In ISchnee fc Goodmaru (120051 ) we report inner and outer powerlaw exponents of —0.8 and 
—1.8, also using just the 450 and 850 fim data and P = 1.5, but using a slightly different 
position as t he center of the nested spher ical shells, so the comparison of powerlaw exponents 
presented in ISchnee fc GoodmanI (120051 ) and those present here would be unfair. Also note 
that the central density of TMC-IC that we derive from the 450, 850 and 1200 /im maps is 
a factor of ~5 lower than when derived using just the 450 and 850 fim maps (and assuming 
f3 = 1.5), which shows the significant adjustments that can result from utilizing additional 
emission maps. An improved value of the central density can make a significant impact on 
the predictions of the dynamical state of a core, and on chemical models. 



6. Summary 



We have used SCUBA data at 450 and 850 fim and MAMBO data at 1200 /xm to 
create maps of the dust temperature an d column density in TMC-IC, improving the results 
presented in lSchnee &: GoodmanI (120051 ). In addition, we are able to estimate the emissivity 
spectral index, finding a value of j3 = 1.8 ± 0.5, based on calibration uncertainties. 

Our analysis shows that noise and calibration errors in maps at 450, 850 and 1200 fim 
would have to be less than ~2% to accurately measure the dust temperature, emissivity 
spectral index and column density from three emission maps. Although such low levels of 
noise and calibration uncertainties are not achievable with the current generation of bolome- 
ters, we show that the dust temperature, emissivity spectral index and column density can 
be accurately mapped if they are fitted at four wavelengths, for instance by including the 
350 yum SHARCII waveband. Obtaining accurate maps of T^^, (3 and Nh are necessary for 
accurate determinations of density and temperature profiles as well as the core mass, which 
determine the time evolution and chemistry of the core. 

The next generation detectors on the JCMT (SCUBA-2) and on APEX (LABOCA), will 
not need to sky-chop, allowing more large-scale structure to be visible, improved calibration 
and fewe r image artifacts, making more accurate determinations of dust properties in cores 
possible JEllislbood : TGiisten et al.lbooeh . 
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Fig. 1. — The true and calculated modified blackbody spectra of material with Ay = 50, 
(3 = 1.5 and Td = 10 K, observed with a 14" beam. The black curve shows the true spectrum, 
and the black crosses show the flux at 450, 850 or 1200 /xm. The blue and red crosses show 
the flux overestimated or underestimated, respectively, by 20%. The blue and red curves 
show the spectrum of dust that passes through the blue/red cross and the two black crosses. 
The dust temperature, emissivity spectral index and column density that would be derived 
from the three data points are also shown. 
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Fig. 2. — The derived emissivity spectral index (/?), the ratio of the derived column density 
to the true column density and the derived temperature as functions of the ratio of the 
measured flux to the true flux at one wavelength (450, 850 or 1200 /xm), with the other 
two fluxes measured without error. The black and red curves show the results for dust with 
Trf = 10 K and = 15 K, respectively. /? = 1.5 in both cases. The curves are cut off where 
the derived temperature is greater than 100 K. 




Fig. 3. — Each panel shows the la spread in the derived column density, temperature or 
emissivity spectral index as a function of the percent noise in all three flux maps at 450, 850 
and 1200 fim. A 5% error in the measured flux corresponds to a signal to noise of 20. The 
black line shows the median value for the derived parameter, and the greyed area shows the 
1 a spread in that parameter for dust at 10 K. The blue (dotted) and red (dashed) hues 
show the median value and 1 a spread for dust at 15 and 20 K, respectively, and each panel 
is repeated for emissivity spectral index (/5 = 1.0, 1.5 or 2.0). Nh, and P are derived as 
in Section [21 The column density is plotted as the ratio of the derived Ay to the input Ay, 
and the temperature and emissivity spectral index are plotted as the difference between the 
derived parameter and the true value. 
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Fig. 4. — Same as Figure [3l except Nfj, and (3 are fit from /crnr observations at 350, 450, 
850 and 1200 /im. 
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Fig. 5. — Each panel shows the 1 a spread in the derived T^, (3 and Nh for dust at = 15 
K and /5 = 1.5 from a least squares fit between the predicted and given fiuxes at 450, 850 
and 1200 /xm. The grey area uses the correct assumption that (3 = 1.5 or Td = 15, while 
the blue (dotted) and red (dashed) lines incorrectly assume that /3 = 1.0 and /5 = 2.0, or 
Td = 10 and = 20, respectively. 
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Fig. 6. — The top row of panels shows the temperature and column density maps of a 
model starless core similar to the central 2' of TMC-IC. The second row shows the predicted 
emission maps at 450, 850, and 1200 fim, using a constant (3 = 1.5. We then add Gaussian 
random noise to the emission maps equal to the measured noise and artifacts in our 450, 
850 and 1200 /xm maps of TMC-IC (80, 9 and 3 mJy, respectively) and show the resultant 
fluxes in the third row. Using the altered fluxes, we attempt to solve for P, and Nh at 
each point in the map, and show the results in the fourth row from the top. When we use 
the value /3 = 1.5 everywhere along with the altered 450, 850 and 1200 fim fluxes to fit 
and Nh, we derive the maps shown in the bottom row. 
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Fig. 7. — The 450, 850 and 1200 jj,m emission maps of TMC-IC, along with the derived /?, 
Trf and A^^^^ (expressed in units of Ay for convenience). The fluxes are in units of mJy per 
14" pixel. The derived parameters are shown within the white 850 fim contour. The (0,0) 
position is at RA=4:41:35.8 DEC=+26:00:42.5 (J2000). Position 'a' is strong at 450 fim but 
not at 850 or 1200 fim. Position 'b' is strong at 850 and 1200 fim but not at 450 fim and 
position 'c' is strong in all three maps. 
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Fig. 8. — The TMC-IC dust color temperature (top) and column density (bottom), derived 
from a fit to the 450, 850 and 1200 /xm emission maps, assuming that the emissivity spectral 
index is constant at /3 = 1.8. 
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Fig. 9. — The density (top) and temperature (bottom) profiles of TMC-IC. Tlie points in red 
are determined by a fit to the 450, 850 and 1200 /im data, assuming [3 = 1.8 and K230 = 0.009 
and t he black points are d e termin ed using only the 450 and 850 /im data assuming (3 = 1.5, 



as m 



Schnee fc Goodman J2005h and /t23o = 0.005, as assumed here. In both cases. 



we 



assume that the central 2 arcminute radius of TMC-IC can be modeled as nested spherical 
shells of uniform density and temperature. 



